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The detailed study of the selected thermodynamic properties of the superconducting phase in the 
molecular hydrogen under the pressure at 428 GPa has been presented. For the increasing value 
of the Coulomb pseudopotential, fi* £ (0.08,0.15), the following results have been obtained: (i) 
the critical temperature decreases from 179 K to 141 K, (ii) the ratio Ri = A (0) /fcsTc differs 
noticeably from the BCS value: _Ri E (4.71, 3.60); (iii) the electron effective mass is large and grows 
slightly together with the temperature ([ m e/ m e] max = 2.2 for T = Tc)- 

PACS numbers: 74.20. Fg, 74.25.Bt, 74.62. Fj 



I. INTRODUCTION 

At low temperature, hydrogen exhibits the nontrivial 
structural behavior under the pressure (p) [lj], Be- 
low 110 GPa the hexagonal-closcd-packed lattice with 
freely rotating molecules is stable (phase I). At higher 
pressures, up to 150 GPa, the broken symmetry phase 
has been observed (phase II). In the pressure range of 
150-350 GPa the so called phase III exists. The exper- 
imental measurements have proved, that the all listed 
phases do not demonstrate the metallic properties. In 
the pressure range from 400 to 500 GPa, the theoretical 
studies predict the existence of the molecular metallic 
phase (the Cmca crystal structure) @, @. Above 
500 GPa, the molecular metallic phase transforms to the 
Cs-IV monatomic phase Q, 0, @. This phase is 
stable at least up to 802 GPa For the extremely high 
value of the pressure (2000 GPa) Maksimov and Savrasov 
have proposed the simple fee structure [1] . 

The molecular and monatomic metallic form of the hy- 
drogen can be the superconductor with the high critical 
temperature (Tc) @. In particular, the calculated values 
of the critical temperature have been presented in Tab. 
HI We notice, that Tc usually has been obtained by us- 
ing the McMillan formula [l5[ , which represents the weak 
coupling limit of the more elaborate Eliashberg approach 
16]. However, in the case of the metallic hydrogen the 
electron-phonon interaction is strong, hence the McMil- 
lan expression is inappropriate. 

For this reason, we have calculated the critical temper- 
ature with the help of the Eliashberg equations. We have 
considered the case p — 428 GPa. Additionally, we have 
studied precisely the properties of the order parameter 
and the electron effective mass. 

In the paper we have taken into consideration the 
Eliashberg set in the mixed representation [17|. This 
approach allows one to obtain the stable solutions on 



TABLE I: The critical temperature for the selected values of 
the pressure; /x* denotes the Coulomb pseudopotential. 
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"The McMillan formula [lj. 

6 The exact solution of the Eliashberg equations |16|| . 
Trobably unstable. 
d The three-band model. 



the real axis, since the analysis does not involve any 
principal-part integrals with singular integrands. 



II. THE ELIASHBERG EQUATIONS 

The Eliashberg equations in the mixed representation 
have been written in the following form pH ]: 
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The Eliashberg set, Eqs. ([T] and [2]) , gives the follow- 
ing solutions on the real axis: the order parameter func- 
tion </> (uj) and the wave function renormalization factor 
Z (uj). The order parameter is defined by the expression: 
A (uj) = (f> (uj) jZ (uj). On the other hand, the imaginary 
axis functions (<f> m = 4> (iuj m ) and Z m = Z (iuj m )) should 
be calculated by using the equations: 
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where the Matsubara energy is given by: 
(tt/P) (2m— 1) and (3 = (^T)' 1 ; ks is the Boltzmann 
constant. 

The pairing kernel for the electron-phonon interaction 
has the form: 
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The Eliashberg function for the molecular hydrogen 
(a 2 F($T)) has been calculated by L. Zhang et al. [lOl ]. 
The maximum phonon energy (f2 roax ) is equal to 508.1 
meV. 

The function fi* (ui m ) = [i*6 (uj c — \uj m \) describes the 
Coulomb repulsion between electrons; //* denotes the 
Coulomb pseudopotcntial, 6 is the Heaviside unit func- 
tion and uj c represents the energy cut-off (uj c = 3f2 max ). 
For the molecular hydrogen we assume: fx* G (0.08,0.15). 

The symbols N (uj) and / (uj) denote the Bose and 
Fermi functions respectively. 

The Eliashberg equations have been solved for 1601 
Matsubara frequencies (M = 800) by using the numeri- 
cal method presented in the paper [ll[ and [l8j|. In the 
considered case, the functions <f> (uj) and Z (uj) are stable 
for T > 11.6 K. 
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FIG. 1: The critical temperature as a function of the Coulomb 
pseudopotential. 
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FIG. 2: The real and imaginary part of the order parameter 
on the real axis for the selected temperatures. The rescaled 
Eliashberg function is also plotted. 



III. RESULTS 

In Fig.Q]we have presented the dependence of the crit- 
ical temperature on the value of the Coulomb pseudopo- 
tential. The exact numerical solutions of the Eliashberg 
equations have been represented by the black circles. The 
dotted and dashed lines represent the calculation of the 
critical temperature by using the McMillan formula and 
the more elaborated Allen-Dynes expression [19( respec- 
tively. It is easy to see, that McMillan formula lowers 
much Tq in the whole range of the Coulomb pseudopo- 
tential's values, whereas the Allen-Dynes expression pre- 
dicts correctly the critical temperature only for very low 
values of fj,*. For this reason we have modified the clas- 
sical Allen-Dynes expression in order to obtain the ana- 
lytical formula, which reproduces the Eliashberg results 
exactly (the solid line in Fig.[l|. In particular, we have 
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FIG. 3: (A) The order parameter on the complex plane for 
the selected values of the temperature. The lines with symbols 
represent the solutions for ui £ (0, f2 ma x), whereas the regular 
lines correspond to the solutions for ui £ (f2 max ,tu c ). (B) The 
temperature dependence of the order parameter. In the both 
cases we assume /i* =0.1. 



fitted the selected parameters in the Allen-Dynes expres- 
sion (Ai and A2) with the help of 250 numerical values 
of Tc (a**)- The final result takes the form: 
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The parameters A, ^/oj2 and oj\ n are equal to 1.2, 207.5 

meV and 141.9 meV, respectively. 

In Fig. [2] we have shown the order parameter on the 

real axis for the selected temperatures, and \i* =0.1; the 

Eliashberg function is also plotted. On the basis of the 

presented results one can state, that both the real and 

imaginary part of the function A (cj) is plainly correlated 

with the shape of the electron-phonon interaction. This 

effect is especially clearly visible for the low values of the 

temperature. The full form of the order parameter on the 

complex plane has been presented in Fig. [3] (A). We have 

stated, that the A (w) values form the distorted spirals 

which shrink with the growth of the temperature. Basing 

on Fig. [3] (A) it is also possible to notice, that the effective 

electron-electron interaction is attractive (Re[A (lS)\ > 0) 

in the range of the frequencies from zero to ~ 0.9f2 max . 
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FIG. 4: The ratio Ri as a function of the Coulomb pseudopo- 
tential. The inset shows the open form of A (0). 
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FIG. 5: The real and imaginary part of the wave function 
renormalization factor on the real axis for the selected tem- 
peratures. The rescaled Eliashberg function is also plotted. 



Taking into consideration the equation: A (T) = 
Re [A (w = A (T))] we have calculated the dependence 
of the order parameter on the temperature (see Fig. [3] 
(B)). Next, the value of the ratio R\ = can be 

obtained, where A (0) denotes the value of the order 
parameter close to the zero temperature and A (0) ~ 
A (T — 11. 6K). In Fig. 2] we have presented the possible 
values of R\ for /i* € (0.08, 0.15). It is easy to see, that in 
the whole range of the considered values of the Coulomb 
pseudopotential, the ratio R\ differs essentially from the 
BCS prediction; [-Ri] BCS = 3.53 20]. Additionally, in the 
inset in Fig. [4] we have shown the dependence of A (0) on 
fi*. 

The second solution of the Eliashberg equations on the 
real axis has been plotted in Fig.[S](/i* = 0.1). The ob- 
tained results prove, that the function Z (u) also clearly 
senses the structure of the electron-phonon interaction. 
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FIG. 6: The wave function renormalization factor on the com- 
plex plane (fi* — 0.1). The lines with symbols represent the 
solutions for lo € (0, £} max ), whereas the regular lines corre- 
spond to the solutions for uj £ (fi max ,tj c ). The inset shows 
the ratio m*/m e as a function of the temperature. 



The shape of the wave function renormalization factor 
on the complex plane has been presented in Fig. [HI We 
see, that in contrast to the order parameter, the function 
Z (lo) weakly depends on the temperature. 

In the framework of the Eliashberg formalism, on 
the basis of the wave function renormalization factor, 
the temperature dependence of the electron effective 
mass (to*) can be calculated. In particular: m*/m e = 
Re [Z (0)] , where the symbol m e denotes the bare electron 
mass. The values of the ratio to* / m e have been presented 
in the Fig's. [6] inset. According to the presented data, it 
is easy to spot, that to* is high in the full range of the con- 
sidered temperatures and [TO*/m e ] max = 2.2 for T = Tc- 
We notice, that at the critical temperature, the electron 



effective mass is independent of jji* 



IV. SUMMARY 

We have calculated the selected thermodynamic prop- 
erties of the superconducting state in the molecular hy- 
drogen (p = 428 GPa). The numerical analysis has been 
conducted in the framework of the one-band Eliashberg 
formalism for the wide range of the Coulomb pseudopo- 
tential's values: /i* e (0.08,0.15). We have proved, that 
the critical temperature is high even for large values of 
A** (Pc] mi „ = 141 K). Next, it has been shown, that the 
superconducting phase is characterized by high value of 
the dimensionless ratio R\, which differs from the BCS 
result; R\ G (4.71,3.60). Finally, we have observed, that 
the electron-phonon interaction strongly enhances the ef- 
fective electron mass for the temperatures from zero to 
Tc- In particular, the maximum of to* is equal to 2.2m e 
for T — Te- 
la the future we will analyze the superconducting state 
in the molecular hydrogen in the framework of the three- 
gap Eliashberg model. It will be done in order to discuss 
the effect of the multi-band anisotropy 
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